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We study the effects of anisotropic pressure on the properties of spheri- 
caUy symmetric, gravitationally bound objects. We consider the full general- 
CN ' relativistic treatment of this problem and obtain exact solutions for various 

forms of the equation of state connecting the radial and tangential pressures. 

> 

\jf^ ' It is shown that pressure anisotropy can have significant effects on the struc- 

' ture and properties of stellar objects. In particular, the maximum value of 

' 2M/R can approach unity (2M /R < 8/9 for isotropic objects) and the surface 

o ■ 

' redshift can be arbitrarily large. 

o 



I. INTRODUCTION 

A common assumption in the study of stellar structure and evolution is that the interior of a 
star can be modeled as a perfect fluid This perfect fluid model necessarily requires that 
the pressure in the interior of a star to be isotropic. This approach has been used extensively 
in the study of polytropes, including white dwarfs, and of compact objects such as neutron 
stars 1^. However, theoretical advances in the last decades indicate that, in many systems, 
deviations from local isotropy in the pressure, in particular at very high densities, may play 
an important role in determining their properties @,|^. 

The physical situations where anisotropic pressure may be relevant are very diverse. By 
anisotropic pressure we mean that the radial component of the pressure, Pr{r), differs from 
the angular components, pe^r) = p^{r) = pt{r). (That pe{r) = p^{r) is a direct consequence 
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of spherical symmetry.) Of course, spherical symmetry demands both to be strictly a func- 
tion of the radial coordinate. A scalar field with non-zero spatial gradient is an example of 
a physical system where the pressure is clearly anisotropic. [This anisotropic character of 
a scalar field occurs already at the level of special relativity, where it is easy to show that 
Pr — Pt = {d(f)/dr)'^]. Boson stars, hypothetical self-gravitating compact objects resulting 
from the coupling of a complex scalar field to gravity, are systems where anisotropic pres- 
sure occurs naturally 0. Similarly, the energy-momentum tensor of both electromagnetic 
and fermionic fields are naturally anisotropic. Isotropy appears as an extra assumption on 
the behaviors of the fields or of the fluid modeling the stellar interior. 

In the interior of neutron stars pions may condense. It has been shown that due to 
the geometry of the n~ modes, anisotropic distributions of pressure could be considered 
to describe a pion condensed phase configuration 0. The existence of solid cores and 
type P superfluidity may also lead to departures from isotropy within the neutron star 
interior Q]. Since we still do not have a detailed microscopic formulation of the possible 
anisotropic stresses emerging in these and other contexts, we take the general approach of 
finding several exact solutions representing different physical situations, modeled by ansatze 
for the anisotropy factor, pt — Pr- As a general rule, we find that the presence of anisotropy 
affects the critical mass for stability, 2M/R, and the surface redshift, Zg- These physical 
consequences of pressure anisotropy are not new. Previous studies have found some exact 



solutions, assuming certain relations for the anisotropy factor [p|-[T^. Our goal here is to 
extend those results, offering a detailed analysis of the changes in the physical properties 
of the stellar objects due to the presence of anisotropy. Hopefully, our results will be of 
importance in the analysis of data from compact objects, as well as in the study of the 
behavior of matter under strong gravitational fields. 

This paper is organized as follows. In the next section, we set up the equations used 
and the assumptions made in our study. We restrict our exact solutions to two classes, 
investigated in sections III and IV respectively. In section III, after reviewing the results of 
Bowers and Liang P], we obtain several new exact solutions for stars of constant density. 
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In section IV, we examine solutions for the case p(r) oc which has been used to model 



ultradense neutron star interiors [0 . We conclude in section V with a brief summary of our 
results and an outlook to future work. In Appendix A, we demonstrate the equivalence of 
the Tolman and Schwarzschild masses for self-gravitating anisotropic spheres. In Appendix 
B we present the general solution for stars featuring an energy density with a constant part 
and a contribution, as discussed in Section IV, in terms of hypergeometric functions. 



II. RELATIVISTIC SELF-GRAVITATING SPHERES 

We consider a static equilibrium distribution of matter which is spherically symmetric. In 
Schwarzschild coordinates the metric can be written as 

ds^ = e'de - e^dr^ - r^dO^ - sin^ Odcf' , (1) 

where all functions depend only on the radial coordinate r. The most general energy- 
momentum tensor compatible with spherical symmetry is 

Tl^ = diag(p, -pr, -pt, -Pt) ■ (2) 

We see that isotropy is not required by spherical symmetry; it is an added assumption. The 
Einstein field equations for this spacetime geometry and matter distribution are 

( - + 4 ) - 4 = ^^Pr ; (3) 
e-(l."-iAV + i(.f +(^^)=8.p,; (4) 

Note that this is a system of 3 equations with 5 unknowns. Consequently, it is necessary to 
specify two equations of state, such as Pr = Pr{p) and pt = Pt{p)- 



3 



It is often useful to transform the above equations into a form where the hydrodynamical 
properties of the system are more evident. For systems with isotropic pressure, this formu- 
lation results in the Tolman-Oppenheimer-Volkov (TOV) equation. The generalized TOV 
equation, including anisotropy, is 

^ = -(p + Pr)^ + -iPt-Pr) , (6) 
ar 2 r 

with 

1 , m(r) + 4:Trr^pr , . 

2^ ~ r(r - 2m) ' ^ 

and 

m(r) = / Anr'^pdr . (8) 



Taking r = i? in the above expression gives us the Schwarzschild mass, M. [This implicitly 
assumes that p = ioi r > R.] 



A more general formula is the Tolman mass formula 

M = j\2T',-Tp{-g)-nHr. (9) 

The equivalence of the Tolman and Schwarzschild masses for systems with anisotropic pres- 
sure is demonstrated in Appendix A. 

In order to solve the above equations we must impose appropriate boundary conditions. 
We require that the solution be regular at the origin. This imposes the condition that 
m(r) — > as r ^ 0. If is finite at the origin then z/' ^ as r ^ 0. The gradient dpr/dr 
will be finite at r = if {pt — Pr) vanishes at least as rapidly as r when r 0. This will be 
the case in all scenarios examined here. 

The radius of the star is determined by the condition Pr{R) =0. It is not necessary for 
Pt{R) to vanish at the surface. But it is reasonable to assume that all physically interesting 
solutions will have Pr,Pt > for r < R. 
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III. EXACT SOLUTIONS FOR p = CONSTANT 



For p = pq = constant, we can write the generalized TOV equation as 

dpr 



.r 2. . ATTr{pr^ + lprPo + '^) 

^ = -(^*-^^0 Y^J^2 



(10) 



Bowers and Liang 1^ solved the generalized TOV equation by considering the following 



equation of state, 



Pt-Pr = C 



r\pr' + IPrPo + ^) 

1 - fvrpor^ 



(11) 



that is, they considered the term related to the anisotropy to be simply proportional to the 
usual hydrodynamic term on the right hand side. The constant parameter C, which we will 
also use, measures the amount of anisotropy. 

They found that the radial pressure is given by 

(1 - 2m/r)« - (1 - 2M/i?)Q 



Pr = Po 



3(1 - 2M/R)Q - (1 - 2m/r) 



where the total mass M = m^R), R is the radius of the system, and Q 



1 _ 3C 

2 47r 



central pressure is given by 



Pc = Po 



1 - (1 - 2M/RY 



(12) 



. The 



(13) 



3(1 - 2M/i?)Q - 1_ 

An equilibrium configuration exists for all values of 2M/R such that Pc is finite. The critical 
model results for that value of 2M/R such that pc becomes infinite. From eq. 0, this occurs 



when the denominator vanishes, that is, when 

(2M/i?)erit 



^3^ 



(14) 



where ^ = 3C/2. We see that for < ^ < 1, (2M/i?)crit can be greater than 8/9, the 
maximum value for an isotropic configuration [|l^]. For a given po and C, the critical mass 
is 

, 3 



crit 



'327rpo' 



1 



1 2 

-)— 

'3^ 



(15) 



Equation |T5| shows that the critical mass is less than the isotropic value when C < 0. When 
C > the critical mass exceeds the isotropic limit. For a given pq, the maximum value of 
the ratio of the critical mass to the isotropic mass approaches 

M,(e = 1)/M, ^ 1.19 , (16) 

where Ma{^ = 1) corresponds to a configuration uniformly filling up to its own Schwarzschild 
radius. This represents a maximum of 19% increase in the stable mass, comparable with 



results for relativistic models of slowly rotating isotropic stars |T^. The increase in stellar 
mass also affects the surface redshift, 

Zs = {l-2M/R)-^ ~1 . (17) 

For the Bowers-Liang solution, the maximum redshift is, 

2c = 3^ - 1. (18) 



For ^ = we recover the well-known isotropic result Zc = 2. The introduction of anisotropy 
removes the upper limit, since as ^ — > 1, Zc can be arbitrarily large. Thus, in principle, a 
modest amount of anisotropy is capable of explaining surface redshifts greater than 2, an 
intriguing possibility. 

For stars with constant energy density po, we considered solutions to two general ansatze, 
which we label cases I and II for convenience. We begin with case I: 
CASE I: The anisotropy factor is written as 

Pt-p^ = Cr^F{pr, Po) il - gTrpo'^M , (19) 

where C measures the amount of anisotropy and the function F{pr,po) includes 6 separate 
cases, 

F = (Pr^ Prpo; PO^; Pr"^ + PrpG\ Pr^ + pG^\PrpG + Po^) • (20) 

For the sake of brevity, we will only discuss 2 of the 6 possible cases. It should be 
straightforward to obtain solutions for the other cases. 
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Case I.l F{pr,po) = po^ 



When F{pr, po) = Po^ the TOV equation becomes 



dpr 
dr 



Pr^ + hrPo + (I - £)Po 



1 - |7rpor2 



(21) 



Thus we find 



rpr 
Jo 



dpr 



rr 

-At: — 
Jr 1 — 



dr 



-.npor' 



(22) 



Pr^ + IprpQ + (I - £)P0^ 

The solutions naturally divide into 3 subcases, according to the value of the anisotropy 
parameter, which we write as A; = C/(27r). 

When k > —1/9 the radial pressure is given by 

(1 - 3k)[{l - 2m/r)4 - (1 - 2M/R)f 



Pr[r) 



Po 



(23) 



(2 + A){1 - 2M/R)f - (2 - A){1 - 2m/r)f 

where A = (1 + 9k)^. Clearly, when A; = we recover the isotropic limit. In figure 1.1, we 
show the radial pressure as a function of the radius parameterized in terms of fractions of the 
critical configuration, (2M/i?)crit- The solutions are normahzed topr{R) — 0. Alternatively, 
to obtain a solution we can fix po and Pc, the central pressure, and compute the value of r 
for which Pr{R) = 0. 




Fig. 1.1: Radial pressure as a function of radius for the case k = 1/2tt, p = const and 
{pt — Pr) oc parameterized in terms of fractions of the critical configuration {2M/ R)crit- 

The central pressure is given by 
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Pc = PO 

Thus, for k > —1/9 we obtain 



il-3k)[l- {1-2M/R)2] 
(2 + A){1 - 2M/R)^ - (2 - A) 



(24) 



2M\ 



1 



crit 



{2 -A) 



For |A;| ^ 1 we find that, correcting shghtly Bowers and Liang 



(25) 



2M\ ^ 8 

~R"/crit " 9 



-Ins) fc + C(A;2) . 



(26) 



We see that a positive anisotropy, k > 0, leads to a violation of the isotropic limit 
(2M//2)crit = 8/9. However, as A; ^ 1/3, (2M/i?)crit 1 and Pr < 0: there is a maxi- 
mum allowed anisotropy in this case. When k < 0, {2M / R)cnt is always less than 8/9. 
The maximum surface redshift for k > —1/9 is 



2 + A 
2- A 



1 . 



(27) 



For k = —1/9, the radial pressure is 



Pr{r) = -po 



1 



1 1^/ l-2m/r 



- 1 



2 -^^(l-2Af/_R,' 

This matches smoothly the result for k > —1/9, as k ^ —1/9. 
When k < —1/9 the radial pressure is given by 



(28) 



Pr[r) 



Po 



B tan 



, 2 \ B (f I- 2m Ir \ ^ 

arctan — H In ■ — 

b) 4:\ \1-2M/r) 



(29) 



with B = {9\k\ — 1)2. Now, the central pressure is 



Po 



B tan (^arctan (^-|^ - j ln(l - 2M/R) ] - 2 



(30) 



and the critical values of 2M/ R are 

f2M 



\ R / crit 



1 — exp 



4 /2 
— arctan — 
B \B 



(31) 



In figure 1.2, we plot the critical values of 2M/R as a function of the anisotropic parameter 
C. 
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The maximum surface redshift for /c < —1/9 is 

Zc = exp 



TT 2 ,2^ 

— — — arctani— ) 
4 B 



1 . 



(32) 



In figure 1.3, we plot the maximum surface redshift as a function of the anisotropic 
parameter C. We find that for ^4 — > 2, i.e., A; — > 1/3, the surface redshift becomes infinite. 
Thus, in this model, positive anisotropies generate arbitrarily large surface redshifts. 




Fig. 1.2 Critical values of 2M/R as a function of the anisotropic parameter C for the 
case p — const and (pt — Pr) oc p^. 




Fig 1.3: Maximum surface redshift as a function of the anisotropic parameter C for the 

case p = const and {pt — Pr) oc p^. 

Case 1.2: F{pr,po) ^ pi 

This solution also separates into 3 subcases, depending on the value of the anisotropy 

parameter, k = C/{2'k). When k > —1/3, the radial pressure is given by 

(1 - 3A;)[(1 - 2m/r)4 - (1 - 2M/i?)4] 



Pr(r) = po 



(2 + A){1- 2M/i?)t - (2 - ^)(1 - 2m/rp 



(33) 
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where A= {1 + Sk)^. In figure 2.1, we show the radial pressure as a function of the radial 
coordinate for fractions of (2M/i?)crit and k — I/tt. 




Fig 2.1: Radial pressure as a function of radius for the case k — 1/n, p — const and 
{Pt — Pr) parameterized in terms of fractions of the critical configuration {2M/R)crit- 



The central pressure is given by 

(l-3A;)[l-(l-2M/i?)# 



Pc = Po 



(2 + A){1 - 2M/Rp -{2 -A) 
The critical configurations for the anisotropy parameter k > —1/3 are given by 



(34) 





[(2-A)l 


\ rC / crit 


[{2 + A)\ 



(35) 



and the corresponding maximum redshift for these values of k are 

^{2 + Ay' 



{2 -A) 

For k — —1/3, the radial pressure is given by 

1 



- 1 



(36) 



Pr{r) = Po 

and the central pressure is 

Pc = Po 
Thus, for A; = —1/3 we have 



2 - Inf ^-^"'A \ 2 



(37) 



2 + ln(l - 2M/i?) 2 



(38) 



\ It / crit 



1-e" 



(39) 
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The solution for A; < —1/3 is 
Po 



Pr{r) 



B tan 



3(1 - k) 

with S = (3|A;| — 1)2. The central pressure is 
Po 



2 \ B / 1- 2m/r \ 
arctan(-J+-ln(^^-^j 



2 , 



(40) 



Pc = 



3(1- A;) 



2 \ B 

B tan ( arctan ( j - j ln(l - 2M/R) 



- 2 



(41) 



Now, 



2M\ 



— I = 1 - exp 

-ft / crit 



4 ,2, n 

— arctan ( — ) 

B ^B' 2 



(42) 



and the corresponding maximum surface redshift is 



Zc — exp 



TT 2 /2 
- --arctan (- 



1 . 



(43) 



Note that the critical values of 2M/R and surface redshifts for cases I.l and 1.2 are identical, 
up to a shift in A; — > A;/3. In figures 2.2 and 2.3 we plot the critical values of 2M/R and the 
maximum surface redshift as a function of the anisotropic parameter C. 




-202468 
c 

Fig 2.2: Critical values of 2M/R as a function of the anisotropic parameter C for the case 
p — const and (pt — Pr) ocp^. 
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2 

Pr 




Fig 2.3: Maximum surface redshift as a function of the anisotropic parameter C for the case 
p = const and {pt — Pr) oc pf.. 



CASE II 

The second class of exact solutions with constant density follows the ansatze, 

Cj^exp{-j^){pr' + lprP0 + lf) 



Pt~Pr 



Po 



1 - fvrpor^ 



(44) 



and 



Pt - Pr 



C^exp{-i){Pr' + lprP0 + ^) 



(45) 



Po i - ^T^por^ 

where tt, > 2 is an integer. The motivation for this choice of anisotropy comes from boson 
stars 1^, where it is found that the anisotropy factor vanishes at the origin and outside the 
star, reaching a maximum somewhere around the approximate radius of the configuration. 
(Boson stars do not have a sharp boundary between the inside and the outside, as the scalar 
field vanishes exponentially for r > R. One may think of it as a diffuse "atmosphere" around 
the denser stellar core.) 

For all cases considered, we found that there are values of C > for which (2M/i?)crit 
can be greater than 8/9 and can approach arbitrarily large values. Let us focus on the 
case where the anisotropy falls exponentially with distance, as in eq. ^ The case for 
the "Gaussian" anisotropy (eq. ^5]) can be solved by following the same procedure. After 
integration we obtain, 
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Pr[r) 



Po 



Z{r) - 1 
3 - Z{r) 



(46) 



where 



(47) 



For a given n the integral can be easily performed and we obtain an expression for the radial 
pressure. Note that for physically acceptable solutions {i.e., with Pr{r) > 0) the function 
Z{r) must satisfy 1 < Z{r) < 3. In figure 3.1 we show the radial pressure as a function 
of radial distance for various fractions of (2M/i?)crit and n — 2. In figure 3.2, we plot the 
anisotropy factor for the same fractions of {2M/R)crit and n — 2. 
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Fig. 3.1: Radial pressure as a function of radius for the case p — const and pt — Pr oc 
(■|)^exp(— parameterized in terms of fractions of the critical configuration {2M/R)crit- 
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Fig. 3.2: Anisotropy factor as a function of the radius for the same case as in fig. 3.1. 

Of more interest is the behavior at the origin. From the above solution, the central 
pressure can be written as. 
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Pc = PO 



where 



Z{0) - 1 
3-Z(0) ' 



Z{0) = (l- SttpoR^) ^^%xp [%'^''-^^exp{-r'/R)dr' 
^ ^ loR^ Jr 



(48) 



The critical configuration is obtained for Z(0) = 3. For example, for n = 2 we obtain, 



(49) 



2M\ 



R 



— 1 exp 



Vcrit ' 9 



(50) 



Note that, as C — > oo, {2M / R)crit 1: there is no maximum positive anisotropy. On the 
other hand, there is a maximum negative anisotropy, beyond which the central pressure be- 
comes negative (^(0) < 1). This is given by \C\ — g(i_2e-i) — 3-12. In the Introduction, 
we noted that boson stars have negative pressure anisotropy, pt — p,. < 0. It is an interesting 
open question if such "exploding" solutions with negative core pressures represent a new 
kind of instability of bosonic stellar configurations. 

In figure 3.3 we plot the critical mass as a function of anisotropy, for several values of n. 
In figure 3.4 we do the same for the surface redshift. For n — 2 the surface redshift is given 

by, 



3 exp 



3 V 



1 . 



(51) 
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Fig. 3.3: Critical values of 2M/R as a function of anisotropy for the case p 
Pt-pr o(. (;|)" exp -{^) with n = 1, 2, 3, 4. 
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Fig. 3.4: Maximum surface redshift as a function of anisotropy for the case p = const and 
Pt - Pr oc (-^)"exp -(-^), with n = 1,2,3,4. 

IV. EXACT SOLUTIONS FOR p ocl/r^ 

We will now consider anisotropic stellar configurations with the following expression for the 
energy density 



where both a and b are constant. The choice of the values for a and b is dictated by the 
physical configuration under consideration. For example, a = 3/7 and 6 = 0, corresponds 
to a relativistic Fermi gas. If we take a = 3/7 and b ^ then we have a relativistic Fermi 
gas core immersed in a constant density background. For large r the constant density term 
dominates (r^ ^ a/36), and can be thought of as modeling a shell surrounding the core. 

In this section we will consider cases where the pressure anisotropy closely resembles the 
behavior of the energy density. Hence we will take 



with c and d constant. The motivation for these ansatze comes from similar approaches 




(52) 




(53) 



used to model the equation of state of ultradense neutron stars [^Q. 
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We have found it convenient to seek solutions for the metric function z/(r) directly, rather 
than solving the generalized TOV equation. We will then use the known functions A(r) and 
z/(r) to find the radial and tangential pressures. From eqs. ^, and ^ we find 



'v" [v'Y' 



'A' 



4 



e'^-u'{- + ^]e-^-(\ + ^]e-^ + \ = M{pt-Pr). (54) 



4 2r 



Introducing a new variable ?/ = e 2 , eq. |5^ becomes 

'A' 1 



-A 



y 



(55) 



Since e = 1 — 2m(r)/r, using eq. B2| we find 



(56) 



In this section we will define the function/^ (x) = 1 — a — bx"^ to simplify our expressions. 
When 6 = 0, we will write Jq = 1 — a. Using the expression for in eq. |5^ and substituting 
for the pressure anisotropy we find 



br"^ — (1 — CL)r'^ y" + (1 ~ a)ry' — [a — c — dr'^)y = 



(57) 



We give the full solution of eq. ^ with a,b,c,d ^ in Appendix B. Below, we will consider 
solutions obtained from choosing specific values for a, b, c and d. 
CASE I: Stars with no crust {b = d = 0) 

We first consider configurations where the energy density is given by 



P 



(58) 



For this density profile, the total mass is M = aR/2 and 



1 - a. 



(59) 



Since for any static spherically-symmetric configuration we expect {2M/ R)crit < 1, we must 
have a < 1. (Also, the metric coefficient grr becomes infinite when a = 1). A density profile 
with this spatial dependence on the radial coordinate was found to be an exact isotropic 
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solution of the TOV equation for the interior of ultra high-density neutron stars by Misner 
and Zapolsky [0. Assuming that the neutron star core can be modeled as a relativistic 
Fermi gas, i.e., Pr = p(r)/3, they found the density to be given by eq. ^ with a = 3/7. 
We note that the Misner-Zapolsky solution cannot be used to construct a complete star, 
since this would require the radius of the star to be infinite. Here, we want to construct 
stars with finite radii and density given by eq. ^ in the context of anisotropic pressure. 
Thus, we impose boundary conditions such that Pr{R) = 0. We also note that Herrera 
investigated anisotropic solutions with similar energy density, in the context of "cracking," 



when perturbations in the fiuid induce anisotropic stresses in the star |13| , p!4| . However, we 
follow a different approach, focusing on the physical properties of static anisotropic solutions. 
With b = d = 0, eq. ^ reduces to an Euler-Cauchy equation, 

(1 - ayy" - (1 - a)ry' + (a - c)y = 0. (60) 

The solutions of this equation divide into three classes, depending on the value of 

(l + c-2a)5 ^ ^ 

<! - 71^ (61) 

(l-a)2 



Case I.l: g is real 

The solution for y is 

with the constants and A_ fixed by boundary conditions. For the case under considera- 
tion here {b = d = 0), the boundary conditions are 

dr R 

Applying the boundary conditions we find 
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2 4g/o 



and A_ — A+{q^—q). 



(64) 



The radial pressure for this case, after substituting the expressions for A^ and is 



R2q _ 



(65) 



(37o' - 1 + 2qll)m + (1 - 37o' + ^ql^y^i _ 

We note that the boundary conditions automatically guarantee that Pr{R) = 0. The radial 
pressure is always greater than zero provided a < 2/3 and > 4c(l — a). Since by definition 
a > 0, the second condition imphes c > 0. Thus, this model does not allow for negative 
anisotropy. Further, since we are considering the case g > 0, we must impose the condition 
1 + c < 2a. Combining the two inequalities for a and c, we obtain, 2a — 1 < c < a^/4(l — a). 
Since we have 0<a<2/3we find that < c < 1/3. We note that for the anisotropic case 
the maximum value of a is 2/3, corresponding to a 33% increase when compared with the 
isotropic case (a = 3/7). In figure 4 we plot the radial pressure, Pr, as a function of the 
radial coordinate r, for a = 3/7 and several values of c. Note that for this choice of a, the 
inequality c < a^/4(l — a) imposes that c < 0.08 for positive pressure solutions. This can 
be seen in the figure. For larger anisotropies, no static self-gravitating stable configuration 
is possible. 



_ c = 0.07 



_ c ^ 0.09 




0.2 0.4 0.6 0.8 1 



Fig. 4: Radial pressure as a function of r for p and {pt — Pr) oc r ^ and q real. 
For r <^ R, we find 



87rr V = 3/o^ - 1 - 2ql^ 



(66) 
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Choosing a = 3/7 we recover, in the hmit c ^ 0, the Misner-Zapolsky sokition [0, with 
Pr = l/(567rr2) = p/3. 



Case 1.2: q = 

The solution for y with g = is 

y(r) = A,(l)+a(^)ta(^) . (67) 

where the constants As and Bg are determined from the boundary conditions. We find 

As = Io and 5, = i^^. (68) 
The radial pressure is given by 

^vrrV = 3/o^ - 1 + (69) 



1 + 

The radial pressure is positive provided a < 2/3. Since we are considering the case q = 0, 
we must require c = 2a — 1 and thus — l<c<l/3. 



Case 1.3: q is imaginary 

The solution for y{r) with q imaginary is 



y = — < As cos 



r 

win ( — 
R 



sin 



r 

win ( — 



with u = \q\. For As and Bs we find 

As = /o, 
The radial pressure is given by 

87rr V = 3/o^ - 1 



1 - 37"2 

and Bs = 

2ula 



2ull 



1 - 3 - 2mJ2 tan (m In (^) 



(70) 



(71) 



(72) 



(1 - Ml) tan (u In (^)) + 2w/o2 

An analysis of eq. |72| shows that > provided that a < 2/3 for any c. Also, since we 
are considering the case q imaginary here, we require that c < 2a — 1. 
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Case II: Including a Crust (a, b,c 0,d = 0) 



We will now derive solutions of the field equations with the density profile given by eq. 



52| . This density profile is essentially a combination of the two profiles (p = const and 
p oc that we have studied so far. We may think of this situation as modelling an 

ultradense core immersed in a background of constant density; at large distances from the 
core (r^ ^ a /3b), the constant density shell dominates the energy density. As in Case I, we 
impose boundary conditions such that Pr{R) = 0. 
For this density profile, 

e-' = m , (73) 

and the equation for j/ = e§ is 

(^{1 — a)r'^ — 6r^^ y" — (1 — a)ry' + ay = 87rr'^y{pt — Pr)- (74) 
The boundary conditions are 

y{R) = im and y'iR)= \~/jll^y (75) 

Here, we have chosen the ansatz for the anisotropic pressure to be the same as the case with 
6 = 0, that is, we do not include a constant contribution to the anisotropy (c? = in eq. 53). 
The general solution for c? 7^ is given in Appendix B. As with the case above with 6 = 0, 
there are three classes of solutions depending on the value of 

(l-2a + c) ^ 
^= (1-a) • 



Case II. 1 q real 

For q real (1 + c > 2a) the solution is 



y 



A, 





'lo + h{r)' 




_Io + h{R)_ 



Iq + h{r) 
_Io + h{R)_ 



+q 



(77) 



Using the boundary conditions as before we find 
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{3li{R) - 1) , h(R) 



+ 



Aqlo 

The radial pressure is given by 

[3/2(r) - 1 - 2M,(r)] + [3/?(r) 



and A- — A+{q — > —q). 



(78) 



Pr 



(t) 




-2'/ /7o+^(,(r)\ 









1 



\r) \, 



(79) 



We note that when 6 = 0, we recover the solution with no crust. Also, if we take a = 
and 6 7^ 0, we have a new class of anisotropic solutions with constant density, and anisotropy 
proportional to r~^. 



Case II.2: g = 

When g = 0, the solution is 



^^'^ " ' -"^'Kr)^ \r) I [/o - h{r)][Io + h{R)] 



with 



As = Ib{R), and Bs 
Here, the radial pressure is given by 



47n 



(80) 



(81) 



8nr% = 3Ii{ry - 1 - 4/;.(r)/oS,e-2 . 



(82) 



Case II. 3: q is imaginary 



When q is imaginary the solution is 



y=[R 



sin(slnF(r)) cos(slnF(r)) 

; — ; „ , , h -De 



'sin(slnF(i?)) cos (slnF(i?)) 



(83) 



with 



As = 2/b(i?)tan(F(i?))tan(2F(i?)) - ■& — ^^^i^ tan(2F), = Ib{R) - A^ , (84) 

sir, 
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and 



F(x) = sln ^ ■ (85) 

(l-a)2 

We intend to investigate numerically the allowed range of the parameters a, 6, c, and ci, 
as well as the stability of these solutions, in a forthcoming publication. 



V. CONCLUSIONS 

We have presented two broad classes of general relativistic exact solutions of spherically 
symmetric stellar configurations exhibiting anisotropic pressure. Our motivation was to 
explore the changes in the general properties of the stars induced by varying amounts of 
anisotropy. In particular, to each of the solutions we have demonstrated that anisotropy 
may indeed change the critical mass and surface redshift of the equilibrium configurations, 
results which we believe are of interest to the astrophysical community and will stimulate 
further investigation. It is an open question if, indeed, anisotropy is revelant for compact 
objects, as, for example, ultradense neutron stars. We have motivated our results based 
on past work on this subject, where isotropic equations of state are at best a reasonable 
hypothesis. Given that we do not have as of yet a complete understanding of the physical 
processes controlled by strong interactions in ultradense matter it is wise to keep an 
open mind to the possibility that anisotropic stresses do occur. Furthermore, hypothetical 
compact objects made of gravitational bound states of bosonic fields, the so-called boson 
stars, are naturally anisotropic |n5|. These objects have been extensively studied in the 



literature, as they represent an interesting new class of compact objects whose stability 
against gravitational collapse comes from a combination of Heisenberg's uncertainty principle 
and model-dependent self-interactions. They may also have some connection with dark 
matter, if perturbations on a fundamental scalar field lead to instabilities which trigger the 
gravitational collapse of overdense regions. 

We have divided our work into two classes of solutions, those with a constant energy 
density and those with an energy density falling as r~^. Within each of these classes we 
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presented several possible cases, which we hope approximate to some extent possible realistic 
objects, including a combined situation where the star's energy density has an ultradense 
interior (p oc r~^) immersed in a shell of constant density. We intend to study the stability of 
these solutions against radial perturbations, as well as provide a detailed numerical analysis 
of the allowed parameter space, in a future work. We also added two appendices, one 
proving the equality between the Tolman and the Schwarzschild mass formulas and the 
other providing the general solution for the p oc + const case in terms of hypergeometric 
functions. 
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Here we will establish the equivalence of the Tolman and Schwarzschild mass formulas for 
static spherically symmetric spacetimes with anisotropic pressure. The general expression 
for the Tolman mass formula is 
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VII. APPENDIX A 




(86) 



Here we have 




and 



T; 



(87) 



Substituting the values of Tq and in equation |8^ we have 




(88) 



Let 
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h= I Anr^ pe'^ dr and h = 4:iTr'^{pr + 2pt)e^dr ; (89) 
Jo Jo 

then, 

Mt = h + h- (90) 

Defining 

rr 

m{r) = / A'Kr'^pdr , (91) 
Jo 

the Schwarzschild mass is given by 

Ms = m{R) . (92) 

Consider the integral Ii. Performing an integration by parts we find 

Ii = me 2 ~ J ) dr. (93) 

The first term evaluates to M5, since m(0) = and, at r = i?, A + z/ = 0. From the field 
equations it follows that 

(A + z/)' 47rr2(p + p,,) 



2 (r - 2m) 

Substituting this expression in equation we find 



(94) 



h = Ms- ""'U —dr. 95 

JO (r — 2mj 

Next we consider the integral l2- Using the generalized TOV equation we can substitute for 
Pt and write 



/, = 4vrr^e^+^ f 3p. + r ^ + i^-Prr' + ^)iP + Pr)\ 
Jo \ dr (r — 2m) J 



(96) 



We will rewrite this as 

h = ha + hb , (97) 

with 
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2/ \ ^+'^ 

/ Anr {3pr)e 2 dr 
Jo 



'2a 



(98) 



and 



'2b 



R 



2 f dpr ^ {Anprr^ + m){p + pr)\ x+. 
Airr r— 1 -, z — ^ ] e ^ dr . 



dr 



Integrating by parts we find 



'2a 



Airr PrC' 



R rR 



(r — 2m) 



'dpr p{X' + u') 



(99) 



Jo \ dr 2 



e 2 dr . 



(100) 



The first term of this equation is equal to zero, since Pr{R) = 0, and substituting for z/' + A' 



from equation |9J we find that 

rR 



dr 



(r — 2m) 



By adding /2a and l2b, we get 



/2= / A.rH'^^^P±me'^dr. 



\ ('^ ~ 2m) 



Finally Ji plus J2 gives 



(101) 



(102) 



Mt = Ms 



(103) 



This demonstrates the equivalence of the Tolman and Schwarzschild mass formulas for static 
spherically symmetric spacetimes with anisotropic pressures. 



VIII. APPENDIX B 



The solution for eq. 57 with a,b,c,d ^ 0, a ^ 1 and boundary conditions given by eq. 



75 



IS 



r \ 1+9 / r \ 



(104) 



Here r = and F is the hypergeometric function: 



F = F{a, A 7, X = 1 + E 

^1 (7)fe k 



{a)k = a{a + l)...{a + k-l) , (105) 
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and 



with 



/3, 7, x) = F'^ ^F{a + 1, /3 + 1, 7 + 1, x) 
ax 7 



F+ = F(a+, 7+, f), F_ = F+(? ^ -g) 



a+ = -[l + 2g-(— ^)] , 
/3+ = -[l + 2g+(— ^)] , 

1 - 2a + c\5 
1 — a 



Using the boundary conditions we find 

^ 1 - 3im - 2qliiR) + 2im{HF-im' 

2h{K){qF+{K) + R[F+{R)]' - RF+{R){ln[F_{K)]}}'' 



-. and A_ = A , 



The radial pressure is given by 

3Iiir) - 1 + 2qli{r) + 27,^(r)(ln[F+(f)])' 



+ 



(3/,^(r) - 1 - 2g4^(r))^(^)-^n^ + 2^^^) 

1 J_ A- F_(r) /r N-2q 
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